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We use perturbative series expansions about a staggered dimerized ground state to compute the 
ground state energy, triplet excitation spectra and spectral weight for a one-dimensional model in 
which each site has an S = | spin Si and a pseudospin T^, representing a doubly degenerate orbital. 
An explicit dimerization is introduced to allow study of the confinement of spinon excitations. The 
elementary triplet represents a bound state of two spinons, and is stable over much of the Brillouine 
zone. A special line is found in the gapped spin-liquid phase, on which the triplet excitation is 
dispersionless. The formation of triplet bound states is also investigated. 
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I. INTRODUCTION 



There has been considerable interest recently in magnetic systems with orbital degeneracy jl],|]j3](|,|H| • Interesting 
^ physics in such systems includes the possibility of orbital ordering and the influence of this on magnetic order, the 
existence and properties of coupled excitations involving both spin and orbital degrees of freedom, and the relation to 
metal-insulator transitions. Materials in which such effects are believed to manifest themselves include the perovskite 
KCuF 3 §, pure and doped manganites Lai__Sr_Mn0 3 §J|, Na 2 Ti 2 Sb 2 and NaV 2 5 f§, LaTi0 3 @, and others. 

Realistic theoretical models of such materials must start from a rather complex Hubbard type model (see e.g. ||) 
with many parameters. To avoid these complications many authors consider an appropriate strong coupling limit 
and electron density so that the charge degrees of freedom are effectively frozen out. The orbital state, for a system 
i i ' with two degenerate orbitals, can then be represented by a pseudospin Tj. Although, in general, the pseudospin 
interactions will not be rotationally invariant, this is often ignored. Thus much recent work has focussed on the "S-T 
r" 1 ' model" 

> ■ _____ 

H = j 1 j^s i - 8, + j 2 J2 T * • ' • ^'_O s < • S ^( T * • w 

(ij) (ij) <*-> 



both in one-dimension |g,|,|,00,|T|,|T|l and in higher dimensions @[D|T(J. 

In this paper we consider the one-dimensional case which can be written in the form 



H = Kj2l(Si ■ S i+1 +j 2 )(T i ■ T l+1 + ji) — jxjz] (2) 



where ji. 2 = Ji^/K. The constant K sets the energy scale and can be set equal to 1. Thus we are left with a 
"■^J ■ 2-parameter Hamiltonian. 

It is known that, at the point j\ = j 2 = 1/4, the Hamiltonian has SU(4) symmetry, and is exactly solvable, by 
the Bethe ansatz method j^n],]^) . At this point the model is gapless with central charge c = 3 and 3 bosonic 
modes. At the point j\ = j_> = 3/4 the model has an exact ground state Jl_} | in which spin and orbital operators 
form dimerized singlets in a staggered pattern. This observation is based on an equivalent realization of the model 
as a 2-leg ladder coupled by a 4-spin interaction 119] . This "matrix product state" is doubly degenerate and gapped. 
The phase diagram of the model thus contains both gapped and gapless regions. The various T = phases and their 
\ boundaries have been mapped out by Itoi et al. (9), using the DMRG technique. 

A number of interesting questions can also be posed regarding the low energy excitations, both in gapless and gapped 
regions. It is by now well known that the elementary excitations in a single spin-half antiferromagnetic chain are 
S = \ spinons p0| . In a system with periodic boundary conditions spinons are formed in pairs, and excitations form a 
two-spinon continuum. However addition of dimerization to a single chain, or coupling of two chains to form a ladder 
leads to confinement of the spinons into bound S = 1 (triplet or magnon) excitations. Spontaneous dimerization, as 
in the Majumdar-Ghosh model pl|, leads to two degenerate ground states. In such a case the elementary excitation 
is a domain wall, separating the two ground states. This excitation also has S = \ and is also refered to as a spinon. 
There has been much recent study of these issues, including the dynamical transition from spinons to triplets/magnons 
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FIG. 1. (a) The dimerized S-T chain, showing the various interaction terms; (b) decoupled plaquettes forming the unper- 
turbed system; (c) one of two exact ground states at the Kolezhuk-Mikeska point. The solid lines represent S = singlet 
pairs. 



In this paper we consider a generalized S-T model in which staggered dimerization is introduced explicitly into 
the S and T chains, as shown in Figure 1(a). There are two reasons for doing this. The first is to provide a natural 
separation of the Hamiltonian into an unperturbed part and a perturbation, allowing the development of "dimer" 
series expansions ]26| . Secondly, by varying the dimerization parameter we are able to move from regions of integer 
spin excitations to S = \ excitations, and study the dynamical transition from triplet to spinon excitations in this 
model. 

We study the ground state energy, energy gap, and excitation spectrum for general values of j%) within the 
gapped phase. At the point j% = 32 = \ the system has an exact dimer ground state, as shown by Kolezhuk and 
Mikeska |lS| ]. We refer to this as the K-M point. The elementary excitations will then be spinons. We would 
expect spontaneous dimerization to exist away from this point, perhaps throughout the gapped phase. However when 
x < 1 one dimerized state becomes the unique ground state and the spinons will become confined into S — 1 triplet 
excitations. At the K-M point the dimerized state is, in fact, an exact ground state for any x. At this point we are 
able to compute substantially longer series for the triplet dispersion and spectral weight. This allows us to study the 
transition from triplet to deconfined spinons as dimerization 5 = 1 — x vanishes. 

We compute dispersion curves for triplet excitations away from the K-M points, specifically along the lines j\ = j 2 , 
and for ji — 1. In each case there is a parameter value where the dispersion curve is flat, and we identify a special line 
in the (j'i,j2) plane having dispersionless excitations. This must mean that along this line excitations are localized 
objects which do not propagate. 



II. DIMER SERIES EXPANSIONS 



In this paper, we study the generalized version of the S-T model, namely the staggered dimerized spin-orbit model, 
using the dimer series expansion approach ]2q| . It is convenient, for this purpose, to view the system as a 2-leg ladder, 
as shown in Figure 1(a). The Hamiltonian for it can be split as 

H = Hq + xV (3) 

where 

Hq = J\ &2i—l ' &2i + J2 ^ T 2 i • T 2 i+i (4) 

i i 

and 

V = Ji ]T S 2J • S 2m + J 2 ^T 2J _! • T 2J + ■ Si+iXT* • T i+1 ) (5) 

i i i 
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FIG. 2. The rescaled ground state per site with x = 1, ji = ji 



Here 6 = 1 — x is the dimerization parameter, and when x — ► 1, the dimerization become zero, and one recovers 
the original Hamiltonian in Eq. fil Our strategy is then to expand the quantities of interest in powers of x and 
to use standard analysis methods |27[ to evaluate these series at x = 1. The unperturbed Hamiltonian Hq consists 
of decoupled plaquettes, as shown in Figure 1(b). The ground state of each plaquette is a product of S and T 
singlets on the respective bonds and the ground state of Ho is then a direct product over plaquettes. This state is 
shown schematically in Figure 1(c), and has energy Eq — — 3( J\ + J 2 )N /16, where N is nnumber of spins (including 
pseudospin). It remains an exact ground state of the full Hamiltonian at the Kolezhuk-Mikeska (K-M) point j\ — 
32 = 3/4. The perturbation V (Eq. ||) includes all interactions between plaquettes. 

To see that the staggered dimer state is the ground state of H at j\ = ji = 3/4 (for any x) we write V as 



V = S 2l • S 2l+1 (J! + KT 2l ■ T 2l+1 ) +J2 T *i-i ■ T M + KS 2l ^ • S 2l ) 



(6) 



Each of the operators (3/4 + T 2 i ■ T^+j), (3/4 + S21-1 • S 2 i) gives zero when acting on this state and hence it is an 
eigenstate of V with eigenvalue zero, and thus an eigenstate of the full Hamiltonian. 



III. RESULTS 



A. General (ji,j2 



We have obtained series to order x 11 for the ground state energy and to order x 9 for the elementary triplet 
excitation energy. The latter corresponds to excitation of either an S or T singlet dimer to an S or T triplet, followed 
by propagation of this triplet along the ladder. We do not include the scries coefficients here but can provide them 
on request. 

The series are then evaluated at the point x = 1 using integrated differential approximants (27) . Figure || shows 
the ground state energy along the line ji — j 2 , and includes the exact SU(4) and K-M values. The series converges 
well, even into the gapless region j\ = j 2 < 1/4, and the energy varies smoothly, consistent with the expected 
Kosterlitz-Thouless transition at the SU(4) point. 

Figure^ show triplet excitation energies versus wavevector q for the case of zero dimerization (x = 1) firstly along 
the line ji = j 2 , and secondly along the line j 2 = 1. In both cases there is a qualitative change in the dispersion 
curves, with the minimum energy at q — n/2 for small j moving to q = for large j. In each case we identify a critical 
ji for which the dispersion curve is virtually flat. In figure ^(a) this occurs at approximately (ji,j 2 ) — (0.725, 0.725), 
slightly lower than the K-M point, whereas in figure ||(b) the dispersion curve is flat at {ji,j 2 ) — (0.65,1). Figure 
H(a) shows the excitation energy at q — n/2 continuously going to zero as ji = j 2 reduces to 0.25, the SU(4) point. 
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FIG. 3. The elementary triplet excitation spectra along the line ji = ji (a) and the line j% = 1 (b) for the system with x = 1. 



From the dispersion curves we can compute the minimum gap A, and the bandwidth W — E(tt/2) — E(0). These 
are shown in Figure ^ for j\ = 22 (and x = 1). We see that the gap opens at j% — ji = 0.25, increases monotonically 
to the K-M point, and remains essentially constant (i.e. A oc K) thereafter. The bandwith W is zero when the 
dispersion is flat. Negative bandwidth by our convention means E(q — ir/2) < E(q = 0), for larger ji, the bandwidth 
W is proportional to J\. 

Itoi et al. have displayed a ground state phase diagram for the model, and designated various phases by symbols 
I- IV. We reproduce this diagram in Figure |[ where we have separated the gapped phase IV into two parts IV and 
VI. The line dividing the regions is where the triplet excitation is dispersionless. 



B- Ui,h) = (3/4,3/4) 

As mentioned above the staggered dimerized singlet state is the exact ground state for any value of the dimcrization 
5 = 1 — x. This feature makes it possible to derive much longer series in x. This feature was also observed along the 
Shastry-Sutherland line in the J\ — J2 — S chain ]25| , ^| ; we do not have a simple physical argument to explain it. We 
have computed series to order x 25 for both the triplet excitation spectrum and the spectral weight. 

The primary focus here is to investigate the transition from triplet to spinon excitations (deconfinement of spinons) 
as x — > 1. In the disconnected dimer limit x — the elementary excitations are localized triplets on one of the spin 
pairs (S2J-1 • S2i) or (T2, • T2i+i). These triplet excitations develop a dispersion for x 7^ but remain well defined 
for all x < 1. At x = 1 the Hamiltonian has full translation symmetry but, as in the well-known Majumdar- Ghosh 
model this symmetry is spontaneously broken leading to two degenerate ground states, one of which remains the 
ground state for x < 1. In this case the elementary excitations are domain walls with S — h interpolating between 
the degenerate ground states. 

Following Singh and Zheng p5[ | we start from the small x limit, and look for breakup of the triplet excitations as 
x — > 1. The triplet dispersion relation is shown in Figure ^| for various values of x. The curves are quite flat for small 
x, and start to develop some dispersion by x ~ 0.9. We expect that the gap at q = behaves as 

A~ A +c(l-xY (7) 

but the series are not regular enough to confirm this or to estimate the exponent v. The lower edge of the two- 
soliton/spinon continuum for x = 1, given by Kolezhuk and Mikeska |Q is shown as a solid line in Figure |. In the 
small q region our triplet dispersion agrees well with this curve, indicating that there are no soliton bound states at 
small q. However for most of the range of q our triplet dispersion is substantially below the two-soliton continuum, 
indicating the formation of soliton bound states. Kolezhuk and Mikeska |18| give a variational ansatz for a triplet 
excitation representing a spinon bound state. Our triplet excitation has signficantly lower energy, and is the favoured 
triplet state. 
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FIG. 4. The minimum gap and bandwidth for elementary triplet excitations for x = 1, ji = ja. 
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FIG. 6. The triplet dispersion relation for system with ji = j% — 3/4 and various x. The solid curve are the lower edge of 
the two-spinon continuum for x — 1, given by Kolezhuk and Mikeska |18|. 



We have also computed series to order x 25 for the spectral weight A(q), and we show results for various x values in 
Figure [?]. For this model the spectral weight satisfies the relation 

A{q)(l - cos(vr - q)) = A(ir - q)(l - cosg). (8) 

Figure ^a) shows a dramatic change in the spectral weight as x — ► 1, approaching the transition to deconfined spinons. 
At x = 1 the spectral weight remain finite for q c < q < it — q c , but vanishes for q < q c or q > 7r — q c . The critical 
wavenumber is estimated to be about q c ~ 7r/6, and this agrees well the triplet dispersion results shown in Figure]^. 
The region over which A(q) remains finite is substantially larger than that for the J\ — J2 — S chain p5| , indicating 
that spinon bound states are more stable in the present model. A simple Dlog Pade analysis of the spectral weight 
A(q) for q < q c or q > n — q c suggests that A(q) vanishes linearly with (1 — x), in agreement with the analysis in ref. 
p4j . This is shown in Figure 0(b) where we show A(q) as a function of x for various q values. 

C. Two- Triplet bound states 

It is also of interest to investigate the possible existence of bound states of the elementary triplet excitations. We 
have done this, using a recently developed series techniques for multiparticle excitation The series for two-triplet 
bound states have been computed up to order x 6 for singlet bound states, and to order x for triplet and quintet states. 
The reason why the singlet series is computed to only 6th order compared to 9th order for the triplet and quintet states 
is that the singlet has the same quantum numbers as the ground state. Thus a much more elaborate orthogonalization 
method is required to implement the cluster expansion for the singlet. Figure ^ shows the excitation spectrum of the 
staggered dimerized S-T chain for parameters (ji, J2) = (0.5,0.75) and dimerization parameter 8 = 1 — x = 0.5. For 
ji < 32 (as here) the S triplets will have lowest energies, and the figure refers to these. The figure show a two-triplet 
continuum, with a number of bound states below the continuum, near q = tt/2 as well as "antibound" states above 
the continuum. The lowest bound state is an S = 2 quintet, with S = singlet and S — 1 triplet bound states at 
higher energy. The behaviour is qualitatively similar to the case of Jj — J2 — S chain p8| , but the quantum numbers 
of the bound states differ. 

To investigate the binding energy as a function of the dimerization parameter we have repeated this calculation 
for different values of x, and show the results in Figure ^, as a binding energy measured relative to the lower edge of 
the continuum. We show results for q = it/2 only, as at this point the binding energy is largest. The binding energy 
vanishes at x = 0, as expected, since the triplet excitations there do not interact. The Qi bound state energy vanishes 
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FIG. 7. (a). The spectral weight A versus momentum q for the system along the line ji = ja — 3/4 for various x (as shown), 
(b) The spectral weight A versus x for the system along the line ji = j2 = 3/4 and various momentum q. 



linearly at small x, whereas the S\ and T\ states vanish quadratically. The binding energies are greatest for x around 
0.5, and vanish as x — » 1. 

IV. DISCUSSION AND CONCLUSIONS 

We have used series expansion methods to investigate properties of a spin model, the dimerized S-T chain, which has 
received much attention recently in connection with orbital ordering in magnetic systems. By including a staggered 
dimerization in the S and T exchange interactions we are able to derive perturbation expansions about a state which 
is the exact dimerized ground state at the Kolezhuk-Mikeska point (ji, J2) = (3/4,3/4), and these expansions are 
regular and analysable throughout the gapless region of the model's phase diagram. 

As well as the ground state energy, we have obtained dispersion curves for elementary triplet excitations. The latter 
show an interesting feature, where the minimum excitation energy changes from q = tt/2 to q = as Ji,j2 increase, 
and the crossover is marked by dispersionless excitations along a particular line in the phase diagram, which we locate. 
The K-M point lies close to this line, but not on it. At the K-M point, using longer series, we find a definite though 
small dispersion. 

In the case of zero dimerization, in the gapped region of phase diagram, it might be supposed that elementary 
excitations consist of dcconfmed spinons (domain walls). However we show that triplet excitations have lower energy 
than that of two spinons, and hence that spinons will bind, except for low momentum. 

We have also used a recently developed series technique to study 2-triplet excitations, and find a variety of bound 
states below and above the continuum. 

After completing this work we learnt of the paper by Martins and Nienhuis They treat exactly the same 

model as discussed here, but for the case K < 0. They find two points at which the model is exactly solvable: one at 
ji = 22 = — l/4, x = 1, the other at ji = j'2 = —1/2, x — 1/2. These lie within phase I (Fig. ||) where the ground 
state is fully ferromagnetic. 
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FIG. 8. The two-particle excitation spectrum of the staggered dimerized spin-orbit chain with j\ = 0.5, ji = 0.75 and x = 0.5. 
Beside the two-particle continuum (gray shaded), there is one singlet bound state (Si), one triplet bound state (Ti), and one 
quintet bound state (Qi) below the continuum, and one singlet antibound state (S2), one triplet antibound state (T2), and one 
quintet antibound state (Q2) above the continuum. The inset enlarges the region near k = tv/2 so we can see Si, Ti and Q2 
below/above the continuum. 
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FIG. 9. The binding energy Eb/K at q = n/2 versus dimerization x for singlet, triplet, and quintet bound states (Si, Ti 
and Qi) of the staggered dimerized spin-orbit chain with j\ — 0.5 and j'2 = 0.75. Several different integrated differential 
approximants to the series are shown. 
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